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Abstract. - We study slip avalanches in disordered materials under an increasing external load 
in the framework of a fiber bundle model. Over-stressed fibers of the model do not break, instead 
they relax in a stick-slip event which may trigger an entire slip avalanche. Slip avalanches are 
characterized by the number slipping fibers, by the slip length, and by the load increment, which 
triggers the avalanche. Our calculations revealed that all three quantities are characterized by 
power law distributions with universal exponents. We show by analytical calculations and com- 
puter simulations that varying the amount of disorder of slip thresholds and the number of allowed 
slips of fibers, the system exhibits a disorder induced phase transition from a phase where only 
small avalanches are formed to another one where a macroscopic slip appeares. 



There is a large variety of non-equilibrium systems 
which exhibit crackling noise, i.e. they respond to a slow 
continuous external driving in the form of bursts of lo- 
, cal events p] . Examples can be mentioned from earth- 
quakes and fracture of disordered materials 2 ,3 ], through 
Barkhausen noise in ferromagnets [3], to martensitic shape 
memory alloys and plastic deformation of solids [516] . Dur- 
ing the last decade experimental and theoretical investi- 
gations revealed that the probability distributions of the 
characteristic quantities of bursts have scale free behavior 
with universal exponents [J[2][6 . 7 . An intense research 
has been initiated to understand the underlying mecha- 
nism of the observed universality. The investigation of 
simple models which grasp the crucial features of systems 
exhibiting crackling noise proved to be essential. Along 
this line, based on the analogy of the plastic deformation 
and fracture of heterogeneous materials, recently, a mi- 
cromechanical model was introduced in Ref. [7] which can 
reproduce the main features of crackling noise in these 
types of systems with only one tuning parameter. 

In the present Letter we study the emergence of crack- 
ling noise in heterogeneous materials which respond to 
an increasing external load by local rearrangements with 
stick-slip mechanism. We consider a fiber bundle model 
[TTJifB] where over-stressed fibers do not break, instead 
they increase their relaxed length in a slip event until they 
can sustain the load. The system is driven by small load 



increments giving rise to the slip of a single fiber which 
may then trigger an entire avalanche of slip events due to 
load redistribution in the bundle. We show by analytic 
calculations and computer simulations that the load in- 
crement triggering the slip bursts, furthermore, the num- 
ber of slipping fibers and the total slip length of the bun- 
dle are all characterized by power law distributions. We 
demonstrate that the amount of disorder and the total 
number of allowed slips play a crucial role in the system: 
a disorder induced phase transition JTJ[6l[Sj|93 i s obtained 
from a low disorder phase where the system snaps with 
macroscopic bursts to the high disorder one where only 
small avalanches pop up. Our model provides an adequate 
description of the micromechanics of disordered systems 
which store hidden length [15] ■ and it can also be consid- 
ered as the fiber bundle analogue of the Burridgc-Knopoff 
model of earthquakes with an infinite range of interaction 

Our model consists of N fibers assembled in parallel. 
Under an increasing external load a the fibers exhibit a 
linearly elastic behavior characterized by the same Young 
modulus E. The important novel element of the model is 
that when the deformation e of a fiber reaches a thresh- 
old value Eth the fiber does not break. Instead, its re- 
laxed length increases until the load reduces to zero on 
the fiber. The mechanism of relaxation is the slip of the 
fiber end, or it can also be interpreted as the unfolding of 
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Fig. 1: (Color online) (a) Constitutive behavior of the bun- 
dle with exponentially distributed quenched failure thresholds. 
o~(e) tends to an asymptotic linear response preceded by a 
longer horizontal plateau when k max increases. The intersec- 
tion of the asymptotic straight lines with the horizontal axis 
indicate the remaining deformation when unloading the sys- 
tem, (b) Breaking probabilities Pk(s) for k m ax = 3. 



subunits of fibers which provide some stored length [15] , 
The slip thresholds e\ h , i = 1, . . . , N are random variables 
with a probability density p(£th) and distribution function 
P{e t h)- After the slip event the fiber gets sticked again so 
that it can support load and can suffer further slips. The 
load kept by fiber i at a deformation e after slipping reads 
as <Ji = E(e — e' l th ) so that no hardening or softening is 
assumed in the system. When a fiber slips again either 
the same slip threshold is retained (quenched disorder) or 
new threshold values can be drawn from the same prob- 
ability distribution p(eth) (annealed disorder). The total 
number of slip events k max a fiber can suffer is a very 
important parameter of the model which can vary in the 
range 1 < k max < +oo. For the load redistribution follow- 
ing slip events we assume an infinite range of interaction, 
i.e. equal load sharing, which is ensured by the condition 
that the strain e is the same for all fibers. In the present 
paper our analysis is restricted to the case of quenched 
disorder so that the load of fiber i in the bundle after 
h < k m ax slips takes the form Ui = E(e — helfj. Further 
details of the model construction can be found in Ref. [15] 
including also the case of annealed slip thresholds. 

Based on the assumption of equal load sharing the con- 
stitutive equation of the parallel bundle can be obtained 
analytically by integrating the load kept by the subsets of 
fibers with different slip indices k [TS] 

cr(e) = Ee [1 - P{Ee)] 

+ X] / p{Eex)E (e - kei) de x 



e/k„ 



p(Eei)E (e - k max Ex) <fei. 



(1) 



Note that the integrals have to be performed over the 
entire loading history of the bundle. For very large de- 
formations e — > co, practically all fibers have suffered 
kmax slips so that Eq. (fT]) can be rewritten as (j(e) ~ 
Ee — k rnax E j^ kmax p(e 1 )e 1 de 1 where the integral pro- 
vides the average value of the slip thresholds (e^). It 
means that the bundle has an asymptotic linear behav- 
ior with the initial value of the Young modulus, however, 
when unloading the system a — > an irreversible perma- 
nent deformation remains whose maximum £™ ax value is 
proportional to the average slip length (£th) an d the num- 
ber of slip events k ma x allowed e™ ax = k m ax(£th)- We 
note that for a finite bundle of N fibers with quenched 
slip thresholds the constitutive equation Eq. ([T]) can be 
written in a discrete form 



AT 

(7(e) =Ee- {E/N)Y,h(e) 
i=i 
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(2) 



where fej(e) denotes the number of slips suffered by fiber 
i up to deformation e. 

In the explicit calculations we use Weibull distributed 
threshold values £th with the probability density function 
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where the parameter A setting the scale of the thresholds 
is fixed to A = 1 in our entire study. The Weibull expo- 
nent 1 < to < +oo is a very important characteristics of 
the system, which controls the amount of disorder in the 
slip thresholds. Increasing the value of m from 1 to infinity 
the probability density Eq. ([3]) varies from the exponential 
distribution to the delta function of zero width. Figure [T] 
illustrates the constitutive curve cr(e) of the model with ex- 
ponentially distributed quenched slip thresholds (to = 1) 
for several different values of k ma x- It can be seen that 
increasing the maximum number of breakings allowed a 
plastic plateau develops, i.e. the final asymptotic linear 
part of the constitutive curve is preceded by a longer and 
longer horizontal plateau. The slope of cr(e) in the asymp- 
totic regime is equal to the Young modulus E — 1 of fibers. 
Note that the simple form of a(e) in Fig. [lja) is the con- 
sequence of the monotonous behavior of the exponential 
distribution, i.e. varying the value of to and k m ax along 
the plastic plateau (7(e) can have a more complex func- 
tional form which will be explored below. Macroscopic 
failure of the system can be captured in the model by as- 
suming that the fibers break after having suffered k m ax 
slips, which has been studied in Ref. [15]. In the present 
paper we focus on the microscopic stick-slip process of the 
fiber bundle with quenched slip thresholds retaining the 
fibers' stiffness after k ma x slips (no breaking). 

Quasi-static stress controlled loading of the fiber bun- 
dle can be performed by incrementing the external load 
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with a small amount Sa just to provoke the slip of a single 
fiber. Since the external load a is kept constant during 
the slip, the load dropped by the slipping fiber must be 
overtaken by the other ones which can give rise to further 
slip events. This way a single slip induced by the load in- 
crement Sa, can trigger an entire avalanche of slips, which 
increases the macroscopic strain e of the system by the 
amount Se. This jerky microscopic dynamics has the con- 
sequence that the deformation of the bundle has a step- 
wise increase under a quasi-statically increasing external 
load a. We characterize the slip avalanches by their size A 
defined as the number of fibers slipping in the avalanche, 
and by the emerging slip length Se which is the increment 
of the strain e of the bundle. All the three quantities, 
the load increment Sa which triggers the avalanche, the 
avalanche size A, and the slip length Se are random vari- 
ables so that the stick-slip process on the micro-scale can 
be characterized by their probability distributions P(Sa) 
P(A), and P(Se), respectively. Note that under strain 
controlled loading no slip avalanches can arise. 

For simple fiber bundles where fibers break irreversibly 
when the local load surpasses their threshold value, it has 
recently been shown pUlfTTj for the case of equal load shar- 
ing that the size distribution of avalanches P(A) can be 
obtained in a closed analytical form as 
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Fig. 2: Constitutive curves a(e) and average number of induced 
slips a(e) for a fixed kmax = 3 varying the value of m (a,b), 
and for a fixed m = 5.0 varying the value of kmax (c, d). Note 
that ml » 1.918. 
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where P denotes the cumulative distribution of the slip 
thresholds. The functional form of P k (e) is presented in 
Fig. [TJ6) for m = 1 and k max = 3. From the above equa- 
tions one can determine the probability density p^ +1 (e) of 
events that a fiber which has suffered k slips until the de- 
e A[a(e)-\na(e)]^ £ ^ u\ formation e was reached, will slip again due to the strain 
increment de 



Here a(e) denotes the average number of fibers which 
break as a consequence of a single fiber failure induced 
by the external load increment at the deformation e. The 
integration over e is carried out up to the critical point 
e c of the system where catastrophic collapse occurs. The 
dominating contribution to the integral is provided by the 
vicinity of the maximum of the exponent of the integrand 
ip(e) = a(e) — lna(e), which is obtained at a = 1. Af- 
ter Taylor expansion of a(e) and of the exponent ip(e) 
about the maximum the asymptotics of the size distribu- 
tion P(A) reduces to the power law form 



P k k + \e) 



1 



1 



P 



1 



< k < k n 



(7) 



P(A) ~ A~ r . 



(5) 



where p is the original probability density of the slip 
thresholds. When the external load is increased by the 
amount Sa at the deformation e to provoke the slip of a 
single fiber i which has already slipped k times, the strain 
increment Se — Se k arising due to load redistribution after 
the slip can be obtained from the constitutive equation of 
finite bundles Eq. ((2|). Keeping the load a fixed during the 
slip, the strain increment reads as Se k — e\ h /N — e/(kN). 
It follows that the average number of fibers a(e) which 
slip as a consequence of a single slip can be determined as 



The exponent t — 5/2 proved to be universal for a broad 
class of disorder distributions where the macroscopic con- 
stitutive curve a{e) of the system has a single quadratic 
maximum [T0l[TT| fT3] . 

In order to understand the dynamics of slip avalanches 
in our model, first the sequence of slipping events has to be 
analyzed. The probability P k (e) that a randomly selected 
fiber in the bundle has suffered exactly k slips up to the 
deformation e can be obtained analytically as 



a(e) = 7V£ 



fe=0 



<5 £ fcPfc +1 ( e ) which leads to the form 



a{e) 



^2 12 



fe=i 



Ee 

k 2? It 



(8) 



It is important to emphasize that the derivative of the 
constitutive equation a{e) of Eq. ([T]) can be expressed in 
terms of a(e) as 



Po(e) 
Pk(e) 



= 1--P(e), 



(6) 



da 
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E 



P 



Ee 



Ee l Ee 



fc=l 



E[l- a(e)] 



(9) 



1 < k < kn 



which show that the constitutive curve a(e) has extrema 
at locations e c where the average number of induced slips 
becomes unity a(e c ) — 1. It also follows from Eq. (j9]) that 
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at the extremal points of a(s) the constitutive curve <r(e) 
has an inflexion point d 2 a/ds 2 = 0. 

Our calculations revealed that varying the amount of 
disorder m and the number of allowed slip events k max the 
statistics of avalanches exhibits a very complex behavior. 
Starting from Eqs. (|1I8I9|) we can determine analytically 
the phase diagram of the system on the m — k max plane 
which classifies all possible functional forms of the con- 
stitutive curves o(e) and of avalanche size distributions 
-P(A). Writing a(e) in the form a(e) = J2k=T a fc( £ )> ^ 
can be seen that each term afe(e) = (e/k 2 )p(Es/k) has a 
single maximum at the strain e k where e k p(Ee k /k) = — 1 
holds. It follows that if a\ has a maximum at e\ with the 
value a\, then the maxima of the other terms afe(e) are 
placed equidistantly as e c k — ke\ with decreasing values 
a k = a\jk. Due to the overlap of the functions dfc(e), the 
consecutive maxima of a(e) do not coincide with that of 
afe(e), however, the equidistant spacing and the decreas- 
ing sequence survive. For the case of Weibull distributions 
the above analysis results in e k — kX and a k = m/(ke), 
where e is the base of natural logarithm. It can be seen 
that for k m ax = 1, when only a single slip is allowed, 
at the critical Weibull exponent m\ = e the constitutive 
curve (7(e) has an inflexion point at the position e\ where 
a(e) has a maximum with the value a{e\) = 1. Similarly, 
for any k max > 1 one can find an m k value of the 
Weibull exponent, where the constitutive curve has an in- 
flexion point with the properties da/de\e k — and 
cPa/de 2 ^ — 0, where at the same time o(£fc max ) = 1 
and da/de\e'j, = hold. 

The phase diagram of the system is presented in Fig. 
02 where the decreasing line represents the m% curve 
which was determined numerically. Note that for k max 
I the critical Weibull exponent is ml = e and m k — > 1 
holds for k max — > +oo. In order to obtain the asymptotics 
of the size distribution of slip avalanches P(A) analyti- 
cally from Eq. Q for parameters along the m k curve, 
the Taylor expansion of a(e) and tp(e) about e k has to 
be continued beyond the first order terms. Following the 
derivations of Ref. |14] the first non-vanishing terms are 
a(e) ~ 1 + C^e - s% m J 2 and tf(e) ~ 1 + C 2 {e - e% m J\ 
which result in a power law asymptotics -P(A) ~ A T 
with the exponent r = 9/4. A similar behavior was found 
in Ref. [14] , where a different physical mechanism led to a 
similar constitutive curve of the system. 

The parameter regime of m and k max below the m k 
curve of Fig. [3] defines the high disorder phase of the model, 
where u(s) is monotonically increasing da/de > and the 
maximum of a(s) is always smaller than 1. Figure [Ha, b) 
illustrates the constitutive behavior a(e) and the average 
number of induced slips a(e) for k max = 3 varying the 
value of to, where the critical disorder parameter is mg ss 
1.918. Since the minimum value of the derivative da/de 
Eq. © is positive in the high disorder phase, the avalanche 
size distribution P(A) behaves as in simple fiber bundles 
when the loading process was stoped at a deformation e m 
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Fig. 3: (Color online) Phase diagram of the system. The de- 
creasing line indicates the 77i£ max curve which separates the 
POP and SNAP regimes. 



before the critical point of macroscopic failure [TU] . Hence, 
for 777, < nT-k max from Eq. ((4]) the size distribution of bursts 
takes the form 

P(A) - A- T e" [a(e "' ) " 1 " lna(£ "' )1A , (10) 

i.e. the power law regime of exponent r = 9/4 is followed 
by an exponential cutoff, where in our case e m is the posi- 
tion of the inflexion point of the constitutive curve. Since 
in the high disorder phase of the model only relatively 
small avalanches pop up away from the phase boundary, 
we call this phase as POP phase, following the terminol- 
ogy of Ref. [S] ■ Figure @|a) presents the size distribution 
of slip avalanches P(A) obtained by computer simulations 
for k max — 7 at different m values in the range m < m?. 
A high quality power law behavior can be observed with a 
diverging cutoff as approaching the critical point m — > 
in agreement with the above derivation. 

In the low disorder regime, above the m k curve of 
Fig. 02 the constitutive curve can have local maxima along 
the plateau regime. It follows from the above derivation 
that at a given value of k max the number of maxima of cr(e) 
is one if the value of m falls in the interval mi < m < e. 
Under stress controlled loading a macroscopic avalanche 
appears resulting in a horizontal jump when the maxi- 
mum of cr(e) is reached. Consequently, this phase of the 
model is called SNAP phase [6]. At very low disorder 
m > e the constitutive curve has a local maximum already 
at kmax = 1 and further maxima occur with decreasing 
height accompanied by a similar oscillating behavior of 
a(e) as k max increases. This feature can be observed in 
FiguresEJc, d), which present the behavior of a(e) and o(e) 
for m = 5 varying the value of k max . In the SNAP phase 
the distribution of avalanche sizes P(A) is determined by 
the first maximum of cr(e) which has a quadratic shape. 
Consequently, similarly to the case of simple fiber bundles, 
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Fig. 4: (Color online) (a) Size distribution of slip avalanches 
P(A) in the POP phase for k ma x = 7 varying the amount of 
disorder m. A high quality power law behaviour is obtained 
with a diverging cutoff as the phase boundary is approached 
m — > m,j ~ 1.572. (b) Comparison of avalanche size distribu- 
tions in the POP and SNAP phases. Excellent agreement is 
obtained with the analytic predictions. 



-P(A) has a power law functional form P(A) ~ A T with- 
out cutoff regime but with an exponent r = 5/2 higher 
than in the POP phase QHIinjIIIl- Burst size distribu- 
tions of the POP and SNAP phases are compared in Fig. 
|4j&), where nice agreement can be observed with the ana- 
lytic predictions. 

In order to understand the transition from the POP to 
the SNAP phase when the amount of disorder m is varied, 
we further analyze the burst size distribution Eq. (fl"0|) . For 
the specific case of k max = 1 we have a(e m ) — a\, where 
a\ = m/e converging to 1 when the amount of disorder is 
decreased in the POP phase m — > m\ — e. After Taylor 
expanding the terms in the exponential function of Eq. 
(fT0|) about mf, we obtain the form 



P(A) ~ A" 



-A/ Ad 



(11) 



where the characteristic burst size has a power law diver- 
gence 



Ao ~ (mi 



(12) 



with the cutoff exponent v = 2. The exact proof is for 
kmax = 1 but our numerical calculations revealed that the 
analytic results Eqs. (Illll2j) hold for all values of k max in 
the POP phase. In order to numerically verify the above 
analytic predictions, we assume that the cutoff avalanche 
size A is proportional to the average size of the largest 
avalanche A ~ (A max ). Figure^a) presents (A max ) ob- 
tained by computer simulations of a bundle of N = 10 7 
fibers with k max = 7 varying the Weibull exponent m in 
a broad range. It can be seen that approaching the phase 
boundary m£ from the POP phase (A ma:E ) diverges, 
i.e. it exhibits a sharp maximum in the finite system. Fig- 
ure [5{6) presents the same data as a function of the dis- 
tance from the critical point m?, where a power law behav- 
ior is evidenced with an exponent v = 2.05±0.05 in a good 
agreement with Eq. ([T2|) . The results imply that varying 
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Fig. 5: The average size of the largest avalanche (A m£t:E } (a) 
and the average largest slip length (Se max ) (c) as a function of 
the Weibull exponent m below the critical point rrij. The same 
quantities are presented in (b, d) as a function of the distance 
from the critical point (7717 — m). 



the amount of threshold disorder the bundle of stick-slip 
fibers undergoes a disorder induced phase transition from 
the high disorder phase where small avalanches pop up to 
the low disorder one where macroscopic avalanches snap 
the system [H3IS]. 

In spite of the complexity of the behavior of avalanche 
sizes, computer simulations revealed a universal functional 
form for both the distribution of the slip length P(8e) and 
for the load increments P(Sa). The total slip length, i.e. 
the strain increment Se occurred during an avalanche of 

+ ej£), where 



*th 



size A reads as Se 
e^ h denotes the slip threshold of fibers taking part in the 
avalanche. The distribution of slip length P(Se) is pre- 
sented in Fig. [SJa) for several values of m. It can be ob- 
served that P(5e) exhibits a universal power law behavior 



P{Se) - 6s- 



(13) 



where the value of the exponent <j> — 2.25 ± 0.05 was ob- 
tained numerically independently of m and k max . Simi- 
larly to the avalanche size A, the characteristic slip length 
defined as the average value of the largest slip length 
(5e max ) is sensitive to the precise shape of a(e) in the 
vicinity of the extremal points. It can be observed in Fig. 
[5Jc) that (Se max ) has a sharp maximum approaching the 
phase boundary from the POP phase and a power law di- 
vergence of the type of Eq. (fT2)) is evidenced in Fig. [5jc?). 
The critical exponent v sa 2 proved to be the same as for 
the characteristic burst size. Figure[6j6) presents that the 
distribution of load increments has a power law decay 

P(5cr) ~Sa- a , (14) 

where the value of the exponent is universal a = 2, it does 
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Fig. 6: (Color online) Probability distribution of slip length 
P(Se) of avalanches (a) and the distribution of load increments 
P(6o~) which trigger the avalanches (b) for different values of the 
Weibull exponent m. Universal power law behavior is obtained 
in both cases. 



not depend neither on m nor on k max . It has to be em- 
phasized that the asymptotics of the distribution P(Sa) is 
determined by the beginning of the loading process where 
large load increments are required to drive the system. 
The reason of universality is that the low stress regime of 
the loaded system is insensitive to the parameters m and 

k 

">max • 

In summary, we studied the statistics of slip avalanches 
in a fiber bundle model where overstressed fibers can re- 
lax in a series of stick slip events. We showed that on the 
macro-scale the stick-slip mechanism leads to plastic be- 
havior with a permanent deformation remaining after the 
load is released. On the micro-scale single slips induced 
by external load increments trigger bursts which give rise 
to a step-wise strain increase. The distribution of load 
increments and of slip length exhibit a universal power 
law behavior with exponents independent of the model's 
parameters. The size distribution of bursts proved to be 
sensitive to the amount of disorder and to the number of 
fibers' slips. Our calculations revealed that at high enough 
disorder only small avalanches pop up, while at low dis- 
order macroscopic avalanches can snap the system. We 
set up the phase diagram of the model and showed that 
the transition between the POP and SNAP phases is 
analogous to disorder induced phase transitions. Besides 
the theoretical interest, our calculations provide insight 
into the statistics of restructurings of systems with hid- 
den length such as biomaterials. Our study was restricted 
to the case of Weibull distributions where the amount of 
disorder can be represented by the exponent m. Gen- 
eralization to other distributions defined over an infinite 
domain is straightforward using the standard deviation as 
a measure of disorder. The value of the critical exponents 
r, v, </>, and a do not have any dependence on the func- 
tional form of the disorder distribution. 



It has been shown in Ref. [B] that the mode of external 
driving has a crutial effect on the critical non-equilibrium 
steady states in slowly driven bistable heterogeneous sys- 
tems with controllable disorder: changing the driving from 
soft to hard a crossover is obtained from the classical 
order-disorder universality class to the quenched Edwards- 
Wilkinson class of SOC type. I our investigations of the 
avalanche statistics only stress controlled loading was con- 
sidered, which corresponds to the perfectly soft driving of 
Ref. [6]. Our phase diagram of Fig. [3] can be considered as 
an extension of the zero softness part of the phase diagram 
of Ref. [B] with the additional degree of freedom of varying 
the number of allowed slip events under an infinite range 
of interaction. It is very interesting to extend our study 
to vary the mode of driving which is currently in progress. 
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